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Surface  Heat  Flux  from  AVHRR  Ice  Surface  Temperature 

OBJECTIVES 

Leads  are  a  key  component  of  pack  ice.  The  long-term  objective  in  this  series  of  studies  of  sea 
ice  has  been  to  observe  and  understand  the  role  of  leads  in  the  surface  heat  balance  of  the  central  Arc¬ 
tic.  The  specific  objectives  have  been  to: 

1.  Learn  how  to  use  satellite  imagery  to  identify  leads  and  determine  their  fractional  area  coverage,  width, 
orientation,  surface  temperature  and  albedo.  Our  principal  data  sources  have  been  AVHRR  and  Landsat 
images.  Much  of  this  effort  required  us  to  learn  how  to  quantitatively  determine  the  surface  tempera¬ 
ture  and  albedo  using  AVHRR  images,  including  all  of  the  required  calibrations  and  corrections. 

2.  Determine  the  heat  flux  from  leads.  The  sensible  heat  flux  from  winter  leads  is  often  opposite  in  sign 
and  one  or  two  orders  of  magnitude  larger  than  what  is  commonly  found  over  thick  pack  ice.  The  flux 
from  leads  often  reverses  the  regional  average  heat  flux  and  as  a  consequence  changes  the  mean  struc¬ 
ture  of  the  lower  atmospheric  boundary  layer  from  stable  to  unstable. 
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3.  Determine  the  surface  temperature  variability.  Much  of  the  surface  temperature  variability  is  in  leads, 
and  we  have  sought  to  determine  the  spatial  scales  of  this  variability  in  different  ice  conditions,  to 
determine  the  proportion  of  the  variability  lost  by  the  1-km  resolution  in  AVHRR  images,  and  to 
account  for  AVHRR’s  coarse  resolution  in  making  estimates  of  the  surface  sensible  heat  flux. 


ACCOMPLISHMENTS 

Our  principal  accomplishments  have  been  published  in  a  series  of  four  papers,  with  additional 
results  published  in  another  four  pt^rs  (see  below  for  a  complete  list).  We  have  two  more  papers  that 
are  nearing  completion. 

1.  Lindsay,  R.  W.  and  D.  A.  Rothrock,  1994;  Arctic  sea  ice  surface  temperature  from  AVHRR,  J.  Climate. 
7, 174-183. 

Abstract; 

The  surface  temperature  of  Arctic  sea  ice  is  estimated  using  the  infrared  channels  of  the 
Advanced  Very  High  Resolution  Radiometer  (AVHRR)  on  satellites  NOAA-10  and 
NOAA-11.  Temperature  statistics  are  analyzed  for  478  cells  measuring  200  km  square 
distributed  over  the  entire  Arctic  Basin  throughout  1989.  The  images  are  cloud  masked 
manually,  and  the  surface  temperature  of  the  cloud-free  area  is  estimated  using  an  algo¬ 
rithm  specific  to  Arctic  pack  ice.  The  rms  error  of  the  estimate  is  thought  to  be  about 
3.2°C,  largely  due  to  uncertainty  in  cloud  masking  and  the  lack  of  knowledge  of  haze  and 
ice  crystal  precipitation  occurrence.  The  mean  temperatures  of  the  cells  range  from  near 
0°C  in  summer  to  below  -45°C  in  winter.  Monthly  averages  range  down  to  -40°C  for  the 
central  Arctic  and  -29°C  for  the  peripheral  seas.  The  monthly-average  standard  deviation 
within  cells  in  the  central  Arctic  is  highest  in  November  (2.2*’C)  and  drops  to  almost  0°C 
in  the  summer.  It  is  commonly  twice  as  large  in  the  peripheral  seas  as  in  the  central  Arc¬ 
tic.  We  also  formulate  a  computation  of  the  sensible  heat  flux  at  the  surface  based  on  the 
spatial  variations  of  the  surface  temperature  estimated  with  AVHRR.  The  contribution  to 
the  heat  flux  by  the  thin  ice  and  leads  that  contribute  to  the  spatial  variability  is  found  to  be 
8  W  m'^  larger  in  the  central  Arctic  than  found  in  earlier  estimates. 

2.  Lindsay,  R.  W.  and  D.  A.  Rothrock,  1994;  Arctic  sea  ice  albedo  from  AVHRR,  J.  Climate,  7,  1737- 
1749. 

Abstract; 

The  seasonal  cycle  of  surface  albedo  of  sea  ice  in  the  Arctic  is  estimated  from  measure¬ 
ments  made  with  the  Advanced  Very  High  Resolution  Radiometer  (AVHRR)  on  the  polar 
orbiting  satellites  NOAA-10  and  NOAA-11.  The  albedos  of  145  2(X)-km-square  cells  are 
analyzed.  The  cells  are  from  March  through  September  1989  and  only  include  those  for 
which  the  sun  is  more  than  10°  above  the  horizon.  Cloud  masking  is  performed  manu^ly. 
Corrections  are  applied  for  instrument  calibration,  nonisotropic  reflection,  atmospheric 
interference,  narrowband  to  broadband  conversion,  and  normalization  to  a  common  solar 
zenith  angle.  The  estimated  albedos  are  relative,  with  the  instrument  gain  is  set  to  give  an 
albedo  of  0.80  for  ice  floes  in  March  and  April.  Tlie  mean  values  for  the  cloud-free  por¬ 
tions  of  individual  cells  range  from  0.18  to  0.91.  Monthly  averages  of  cells  in  the  central 
Arctic  range  from  0.76  in  April  to  0.47  in  August.  The  monthly  averages  of  the  within- 
cell  standard  deviations  in  the  central  Arctic  are  0.04  in  April  and  0.06  in  September.  The 
surface  albedo  and  surface  temperature  are  correlated  most  strongly  in  March  (R  =  -0.77) 
with  little  correlation  in  the  summer.  The  monthly  average  lead  fraction  is  determined 
from  the  mean  potential  open  water,  a  scaled  representation  of  the  temperature  or  albedo 
between  0.0  (for  ice)  and  1.0  (for  water);  in  the  central  Arctic  it  rises  from  an  average 
0.025  in  the  spring  to  0.06  in  September.  Sparse  data  on  aerosols,  ozone,  and  water  vapor 


2 


in  the  atmospheric  column  contribute  uncertainties  to  instantaneous,  area-average  albedos 
of  0. 1 3, 0.04,  and  0.08.  Uncertainties  in  monthly  average  albedos  are  not  this  large.  Con¬ 
temporaneous  estimation  of  these  variables  could  reduce  the  uncertainty  in  the  estimated 
albedo  condsiderably.  The  poor  calibration  of  AVHRR  channels  1  and  2  is  another  large 
impediment  to  making  accurate  albedo  estimates. 

3.  Lindsay,  R.  W.  and  D.  A.  Rothrock,  1995:  Arctic  sea  ice  leads  from  Advanced  Very  High  Resolution 
Radiometer  images,  J.  Geophys.  Res.,  100  (C3),  4533-4544. 

Abstract: 

A  large  number  of  advanced  very  high  resolution  radiometer  (AVHRR)  images  from 
throughout  1989  are  analyzed  to  determine  lead  characteristics.  ITie  units  of  analysis  are 
square  2(X)-km  cells,  and  there  are  270  such  cells  in  the  data  set.  Clouds  are  masked  man¬ 
ually.  Leads  are  determined  from  images  of  the  potential  open  water  5,  a  ^aled  version  of 
the  surface  temperature  or  albedo  that  weights  tWn  ice  by  its  thermal  or  brightness  impact. 

The  lead  fraction  is  determined  as  the  mean  5;  the  monthly  mean  lead  fraction  rmges  from 
0.02  in  winter  to  0.06  in  summer  in  the  ceiitral  Arctic  and  is  near  0.08  in  the  winter  in  the 
peripheral  seas.  A  method  of  accounting  for  lead  width  sampling  orors  due  to  the  finite 
sample  areas  is  introduced.  In  the  central  Arctic  the  observed  mean  lead  width  for  a 
threshold  of  5  =  0. 1  ranges  from  2  or  3  km  (near  the  resolution  of  the  instruinent)  in  the 
winter  to  6  km  in  the  summer.  In  the  peripheral  seas  it  is  about  5  km  in  the  winter.  Width 
distributions  are  often  more  heavily  weighted  in  the  tail  than  exponential  distributions  and 
are  well  tqjproximated  by  a  power  law.  The  along-track,  number  density  power  law  N  - 
aw*  has  a  mean  exponent  of  fc  =  1.60  (standard  deviation  0.18)  and  shows  some  seasonal 
variability.  Mean  floe  widths  in  the  central  Arctic  are  40  to  50  km  in  the  winter,  dropping 
to  about  10  km  in  the  summer.  For  floes  the  power  law  has  a  mean  exponent  of  0.93  and 
exhibits  a  clearer  annual  cycle.  Lead  orientation  is  determined  with  a  method  based  on  the 
direction  of  maximum  extent 


4.  Lindsay,  R.  W.,  D.  B.  Percival,  and  D.  A.  Rothrock,  1996:  The  Discrete  Wavelet  Transform  and  the 
Scale  Analysis  of  the  Surface  Properties  of  Sea  Ice,  IEEE  Trans.  Geoscience  and  Remote  Sensing,  in 
press. 


Abstract: 

The  formalism  of  the  one-dimensional  discrete  wavelet  transform  (DWT)  based  on 
Daubechies  wavelet  filters  is  outlined  in  terms  of  finite  vectors  and  matrices.  Both  the 
scale-dependent  wavelet  variance  and  wavelet  covariance  are  considered  and  confidence 
intervals  for  each  are  determined.  The  variance  estimates  are  more  accurately  determined 
with  a  maximal-overlap  version  of  the  wavelet  transform.  The  properties  of  several 
Daubechies  wavelet  filters  and  the  associated  basis  vectors  are  discussed.  BoA  the  Mallat 
orthogonal-pyramid  algorithm  for  determining  the  DWT  and  a  pyramid  algorithm  for 
determining  the  maximal-overlap  version  of  the  transform  are  presented  in  terms  of  finite 
vectors.  As  an  example,  we  investigate  the  scales  of  variability  of  the  surface  temperature 
and  albedo  of  spring  pack  ice  in  the  Beaufort  Sea.  The  data  analyzed  are  from  individual 
lines  of  a  Landsat  1^  image  (25-m  sample  interval)  and  include  both  reflective  (chaimel 
3,  30-m  resolution)  and  thermal  (channel  6, 120-m  resolution)  data.  The  wavelet  variance 
and  covariance  estimates  are  presented  and  more  than  half  of  the  variance  is  accounted  for 
by  scales  of  less  than  800  m.  A  wavelet-based  technique  for  enhancing  the  lower-resolu¬ 
tion  thermal  data  using  the  reflected  data  is  introduced.  The  simulated  effects  of  poor 
instrument  resolution  on  the  estimated  lead  number  density  and  the  mean  lead  width  are 
investigated  using  a  wavelet-based  smooth  of  the  observations. 
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The  two  p^)a'S  that  are  nearing  completion  are: 


5.  Rothrock,  D.  A.  and  R.  W.  Lindsay:  Regional  average  heat  flux  from  leads 
Abstract 

The  regionally  averaged  turbulent  and  radiative  heat  fluxes  from  leads  are  determined 
using  high-resolution  satellite  and  aircraft  observations  of  the  surface  temperature.  We 
determine  the  heat  flux  from  leads  by  extracting  sets  of  one  dimensional  transects  from 
two  dimensional  thermal  images  of  sea  ice  or  by  using  surface  temperature  measurements 
from  aircraft  transects.  The  flux  from  the  leads  is  found  by  comparing  the  lead  surface 
temperature  to  the  tempCTature  of  the  nearby  floes.  A  fixed  wind  spe^  is  assumed.  We 
outline  the  details  of  the  procedures  used  for  determining  the  sensible,  latent,  and  radiative 
heat  fluxes  from  leads  and  the  propaties  of  the  heat  transfer  coefficients  used  in  these  cal¬ 
culations,  we  then  use  this  proc^ure  on  eighteen  Landsat  scenes  and  two  aircraft  transects 
and  then  determine  the  spatial  scales  important  for  lead  heat  fluxes.  The  spatial  scales  that 
are  important  for  the  flux  from  leads  are  determined  using  low-pass  filtared  versions  of  the 
data  with  filto'S  based  on  the  orthogonal  discrete  wavelet  transform.  A  correction  proce¬ 
dure  for  low-resolution  data  is  determined.  Finally,  we  estimate  the  heat  flux  from  leads 
for  114  AVHRR  scenes  from  throughout  the  cold  seasons  of  a  single  year. 


6.  Lindsay,  R.  W.:  Spatial  variability  of  the  surface  temperature  of  sea  ice. 

Abstract 

The  spatial  scales  of  the  variability  of  the  surface  temperature  of  sea  ice  are  determined 
with  the  use  of  the  discrete  orthogonal  wavelet  transform.  The  wavelet  variance  is  deter¬ 
mined  for  18  Landsat  TM  images  from  both  the  Arctic  and  the  Antarctic  for  scales  ranging 
from  30  m  to  30  km.  Aircraft  and  AVHRR  data  are  also  analyzed,  extending  the  range  of 
scales  observed  to  1  m  to  100  km.  The  degree  of  spatial  concentration  of  the  variance  is 
determined  and  the  relationships  between  the  sensible  heat  flux  and  the  spatial  variability 
are  determined.  The  lead  width  distribution  and  the  lead  width/lead  temperature  bivariate 
distributions  are  determined  and  compared  to  various  parameterized  distributions. 
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Papers  in  progress  2 
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